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In this paper we define a new class of weighted complex networks sharing several properties with 
fractal sets, and whose topology can be completely analytically characterized in terms of the involved 
parameters and of the fractal dimension. The proposed framework defines an unifying general theory 
of fractal networks able to unravel some hidden mechanisms responsible for the emergence of fractal 
structures in Nature. 

PACS numbers: 64.60.aq Complex Networks, 89.75.Fb Structures and organization in complex systems, 
89. 75. Da Scale-free networks, 05.45.Df Fractals 

I. INTRODUCTION 

Complex networks have recently attracted a growing interest of scientists from different fields of research, mainly 
because complex networks define a powerful framework for describing, analyzing and modeling real systems that 
can be found in Nature and/or society. This framework allows to conjugate the micro to the macro abstraction 
levels: nodes can be endowed with local dynamical rules, while the whole network can be though to be composed by 
hierarchies of clusters of nodes, that thus exhibits aggregated behavior. 

The birth of graph theory is usually attributed to L. Euler with his seminal paper concerning the "Konigsberg 
bridge problem" (1736), but it is only in the 50's that network theory started to develop autonomously with the 
pioneering works of Erdos and Rcnyi [ij . Nowaday network theory defines a research field in its own 0, II] and the 
scientific activity is mainly devoted to construct and characterize complex networks exhibiting some of the remarkable 
properties of real networks, scale-free0, small-world communities Q, just to mention few of them. 

In a series of recent papers 0, H, [^authors proposed a new point of view by constructing networks exhibiting 
scale-free structures following ideas taken from fractal construction, e.g. Koch curve or Sierpinski gasket. The aim of 
the present paper is to generalize these latter constructions and to define a unifying theory, hereby named Weighted 
Fractal Networks, WFN for short, whose networks share with fractal sets several interesting properties, for instance 
the self-similarity. 

The WFN are constructed via an explicit algorithm and we are able to completely analytically characterize their 
topology as a function of the parameters involved in the construction. We are thus able to prove that WFN exhibit 
the "small-world" property, i.e. slow (logarithmic) increase of the average shortest path with the network size, and 
large average clustering coefficient. Moreover the probability distribution of node strength follows a power law whose 
exponent is the Hausdorff (fractal) dimension of the "underlying" fractal, hence the WFN are scale-free. 

WFN also represent an explicitely computable model for the renormalization procedure recently applied to complex 
networks 

The paper is organized as follows. In the next section we will introduce the model and we outline the similarities 
with fractal sets. In Section |lTT] we present the analytical characterization of such networks also supported by dedicated 
numerical simulations. We then introduce in Section HVl a straightforward generalization of the previous theory, and 
thus we conclude by showing a possible application of WFN to the study of fractal structures emerging in Nature. 

II. THE MODEL 

According to Mandelbrot [l^ "a fractal is by definition a set for which the Hausdorff dimension strictly exceeds the 
topological dimension" . One of the most amazing and interesting feature of fractals is their self-similarity, namely 
looking at all scales we can find conformal copies of the whole set. Starting from this property one can provide rules 
to build up fractals as fixed point of Iterated Function Systems [T^ . [Tsj . IFS for short, whose Hausdorff dimension 
is completely characterized by two main parameters, the number of copies s > 1 and the scaling factor < / < 1 
of the IFS. Let us observe that in this case this dimension coincides with the so called similarity dimension p^ . 

dfract = - logs/ log/. 
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G' =Ts.f.a G 



FIG. 1: The definition of the map Tsj^a- On the left a generic initial graph G with its attaching node a (red on-line) and a 
generic weighted edge w £ G (blue on-line). On the right the new graph G' obtained as follows: Let G'^', . . . , G'°' be s copies 
of G, whose weighted edges (blue on-line) have been scaled by a factor /. For i = 1, . . . , s let us denote by a*-*' the node in 
G^'' image of the labeled node a G G, then link all those labeled nodes to a new node a' (red on-line) through edges of unitary 
weight. The connected network obtained linking the s copies G'*' to the node a' will be by definition the image of G through 
the map: G' = %j,a{G). 



The main goal of this paper is to generalize such ideas to networks, aimed at constructing weighted complex 
networks [20| with some a priori prescribed topology, that will be described in terms of node strength distribution, 
average (weighted) shortest path and average (weighted) clustering coefficient, depending on the two main parameters: 
the number of copies and the scaling factor [2l|. Moreover taking advantage of the similarity with the IFS fractals, 
some topological properties of the networks will depend on the fractal dimension of the IFS fractal. 

Let us fix a positive real number / < 1 and a positive integer s > 1 and let us consider a (possibly) weighted 
network G composed by N nodes, one of which has been labeled attaching node and denoted by a. We then define a 
map, Tsj_a, depending on the two parameters s, / and on the labeled node a, whose action on networks is described 
in Fig. [i 

So starting with a given initial network Go we can construct a family of weighted networks {Gk)k>o iteratively 
applying the previously defined map: Gk 7^./,a(G/c_i). 

Because of its general definition, the map Ts.f\a improves the constructions recently proposed in 0, allowing 
us to consider all possible IFS fractals in a unified scheme instead of using "ad hoc" constructions. For the sake of 
completeness we present numerical results for two WFN: the Sierpinski one (see Fig. [5]) and the Cantor dust (see 
Fig. [31). 

Given Go and the map 7^,/, a we are able to completely characterize the topology of each Gk and also of the limit 
network Goo, defined as the fixed point of the map: Goo = Xi.f.aiGoo). Thus the WFN undergo through a growth 
process strictly related to the inverse of the renormalization procedure \idl . [ll| ; at the same time Goo will be infinitely 
renormalizable . 



III. RESULTS 



The aim of this section is to characterize the topology of the graphs Gk for all fc > 1 and Goo, by analytically 
studying their properties such as the average degree, the node strength distribution, the average (weighted) shortest 
path and the average (weighted) clustering coefficient. 

At each iteration step the graph Gk grows as the number of its nodes increases according to 

Nk = sf^No + (.s^ - l)/(s - 1) , (1) 

being Nq the number of nodes in the initial graph, while the number of edges satisfies 

Ek = s'^Eo + s(s^ - l)/(s - 1) , (2) 

being Eq the number of edges in the graph Gq. Hence in the limit of large k the average degree is finite and it is 
asymptotically given by 

Efe ^ s + Eojs - 1) 

Nk 1 + (s - l)iVo ■ ^ ' 
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FIG. 2: The "Sierpinski" WFN, s = 3, / = 1/2 and Go is composed by a single node. From the left to the right Gi, G2, G3 
and G4. Gray scale (color on-line) reproduces edges weights: the darker the color the larger the weight. The dimension of the 
fractal is log3/log2 ^ 1.5850. Visualization was done using Himmeli software [l7[. 
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FIG. 3: The "Cantor dust" WFN, s = 4, / = 1/5 and Go is a triangle. From the left to the right Go, Gi, G2 and G3. 
Gray scale (color on-line) reproduces edges weights: the darker the color the larger the weight. The dimension of the fractal is 
log 4/ log 5 ~ 0.8614. Visualization was done using Himmeli software 



Let us denote the weighted degree of node i € Gk, also called node strength [l^ , by wp'' = J2j '^ij'' j being w^'^'' 
the weight of the edge (ij) e Gk', then using the recursive construction, we can explicitly compute the total node 
strength, Wk = ^'^d, provided sf ^ 1, easily show that 

Because / < 1, we trivially find that the average node strength goes to zero as k increases: Wk/ N}^ — > 0. 

k — ^00 
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FIG. 4: Node Strengths Distribution. Plot of the renormalized node strengths distribution logj^Q [gk{x)/ D], where D = dfrac in 
the homogeneous case, while D = ~ iog(/°^"/ ) non-homogeneous one. Symbols refer to : □ the finite approximation G14 

with 2391484 nodes of the "Sierpinski" WFN, s = 3, / = 1/2 and Go is formed by one initial node; Q the finite approximation 
Gil composed by 3495253 nodes of the "Cantor dust" WFN, s = 4, / = 1/5 and Go is made by a triangle; A the finite 
approximation Gn composed by 3495253 nodes of the non-homogeneous "Cantor dust" WFN, s = 4, fi = 1/2, (2 = 1/3, 
/a = 1/5, /4 = 1/7 and Go is formed by a triangle. The reference line has slope —1; linear best fits (data not shown) provide a 
slope -0.9964 ± 0.034 with = 0.9993 for the Sierpinski WFN, a slope -1.002 ± 0.064 with = 0.9996 for the Cantor dust 
WFN and a slope -1.006 ± 0.024 with R'^ = 0.9976 for the non-homogeneous Cantor dust WFN. 



A. Node strength distribution. 



Let gk{x) denote the number of nodes in Gk that have strength ujf^^ = x and let us assume go to have values in 
some finite discrete subset of the positive reals, namely: 

g(){x) > if and only if x G {xi, . . . ,Xm \ , 

otherwise go{x) ~ 0. Using the property of the map Tsj.a we straightforwardly get guix) = sgk-i{x/ f) provided [l^l 
X ^ s and x ^ fs + I, from which we can conclude that for all k: 

gk{x) = s^go{x/f), gk{fs + l)^s and gk{s) ^ I . (4) 

This implies than the node strengths are distributed according to a power law with exponent djract = "logs/ log/, 
that equals the fractal dimension of the fractal obtained as fixed point of the IFS with the same parameters s and /. 
In fact defining Xik — f^Xi we get: 

log gk{xik) = /clogs + log 50 

logs , , , .X logs 

= 7 log^ifc + log5o(a;j) - 1 jlogxi, 

log/ log/ 

namely (see Fig. U]) 

gk{x) ^ C/x"'- . (5) 



B. Average weighted shortest path. 



By definition the average weighted shortest path [J of the graph Gk is given by 
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FIG. 5: The average weighted shortest path. Plot of the renormaUzed average weighted shortest path Afe versus the iteration 



number k, where Afe 



(s-F)(s^-F) 
2s^{s-l) 



and F = fi + ■ ■ ■ + fs for the non-homogeneous case, while F = sf for the homogeneous 
3, / = 1/2 and Go is formed by one initial node; Q the "Cantor dust" 

4, /i = 1/2, /2 = 1/3, 



one. Symbols refer to ; □ the "Sierpinski" WFN, s 

WFN, s = 4, / = 1/5 and Go is made by a triangle; A the non-homogeneous "Cantor dust" WFN, s 
fs = 1/5, /4 = 1/7 and Go is formed by a triangle. 



where 

ijeGk 

(k) 

being plj ' the weighted shortest path linking nodes i and j in Gk ■ 

To simplify the remaining part of the proof it is useful to introduce A^°'°'' = X^ieG^ ^'ial' ^'^^ ^^'■^ '^'^ weighted 
shortest paths ending at the attaching node, Ofc £ Gk- One can prove (see Appendix lA ip that for large k the 
asymptotic behavior of A^"'''' is given by 

.(a.) iVo(s-l) + l 

Using the construction algorithm and its symmetry one can prove (see again the Appendix lA 2\ that A^. satisfies 
the recursive relation 

Afc = sfKk-i + 2s[(s - l)Nk-i + l][Nu-i + fk^CX^ > (9) 
that provides the following asymptotic behavior in the limit of large k (see Fig. [5]) 

A = ^ 2(.-l) 

Arfc(A^fc-l) fe^^^ (l-/)(s-/)- ^ ' 

We can also compute the average shortest path, £k, formally obtained by setting / = 1 in the previous formulas ([B]) 
and ([7]) . Hence slightly modifying the results previously presented we can prove that asymptotically we have 

4 - 2(k-^)^ ^logNk, (11) 

fc^oo Y S — 1 J k^oo logs 

where the last relation has been obtained using the growth law of Nk given by equation ([T|) (see Fig.[S]). Let us remark 
that the average shortest path is a topological quantity and thus it doesn't depend on the scaling factor, that's why 
we don't report in Fig. [H] the case of the non- homogeneous WFN. 

Thus, as previously stated, the network grows unbounded but with the logarithm of the network size, while the 
weighted shortest distances stay bounded. 
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FIG. 6: The average shortest path Ik as a function of the network size (semilog graph). Plot of the renormalized average 
shortest patli Ilk versus the network size Nk, where Ik — £k^-^^- Symbols refer to : □ the "Sierpinski" WFN, s — 3, f — 1/2 
and Go is formed by one initial node; Q the "Cantor dust" WFN, s — 4, f = 1/5 and Go is made by a triangle. The reference 
line has slope 1. Linear best fits (data not shown) provides a slope 0.9942 ± 0.019 and R'' = 1 for the Sierpinski WFN and a 
slope 0.9952 ± 0.019 and = 1 for the Cantor dust WFN. 



C. Average clustering coefficient. 

The average clustering coeffieient Q of the graph Gk is defined as the average over the whole set of nodes of 

the local clustering coefficient c^*'', namely < Ck >= Ck/Nk, where Ck = X^ieG^ '^i'^^- Because of the construction 
algorithm the number of possible triangles, hence the local clustering coefficient, at each step increases by a factor s; 
thus after k iterations we will have Ck ~ s'^Co, being Co the sum of local clustering coefficients in the initial graph. 
We can thus conclude that the clustering coefficient of the graph is asymptotically given by: 

s - 1 < Co > TVq 

'^'"'^kZ^ s (s_l)iVo + l- ^ ' 

On the other hand, one can use edges' values to weigh the clustering coefficient [l^; hence generalizing the previous 
relation, we can easily prove that the average weighted clustering coefficient of the graph is asymptotically given by: 

k^oo fs (s-l)iVo + l fe^oo ^l/rf/..c* ' ^ 

where once again, the fractal dimension dfract of the IFS fractal play a relevant role. 



IV. NON-HOMOGENEOUS WEIGHTED FRACTAL NETWORKS 



The aim of this section is to slightly generalize the previous construction to the case of non-homogeneous scaling 
factors for each subnetwork G^*-*. So given an integer s > 1 and s real numbers /i, . . . , /s S (0, 1), we modify the map 
%.f.a by allowing a different scaling for each edge weight according to which subgraph it belongs to: if the edge w^^\ 
image of w S G, belongs to G^^\ then w^^"^ = fjW. 

Let us remark that the construction presented in the former Section |TT] is a particular case of the latter once we 
take /i = • • ■ = /s = /; we nevertheless decided for a sake of clarity, to present it before, because the computations 
involved in this latter general construction could have hidden the simplicity of the underlying idea. We hereby present 
some results for the non-homogeneous "Cantor dust" WFN (see Fig.[7|). 

Using the recursiveness of the algorithm we can, once again, completely characterize the topology of the non- 
homogeneous WFN, moreover only the weighted quantities will vary with respect to the homogeneous case. For 
instance, a straightforward, but cumbersome, generalization of the computations presented in the previous Sections 
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FIG. 7: The non-homogeneous "Cantor dust" WFN, s = 4, /i = 1/2, /2 = 1/3, /a = 1/5, /4 = 1/7 and Go is formed by a 
triangle. From the left to the right Go, Gi, G2 and G3. Gray scale (color on-line) reproduces edges weights: the darker the 
color the larger the weight. Visuahzation was done using Himmeli software [l7j . 



allows us to prove that the average weighted shortest path exhibits the following asymptotic behavior (see Fig. ^ 

2s2(s - 1) 
, — ^ 7 EnTl En ' 

k^oo (s — )(s^ ~ j 

where F = fi + ■ ■ ■ + fg. Let us observe that Eq. reduces to Eq. (fTn|) once we set /i = ••• = /, = / and thus 

Let 50(2^) denote the number of nodes with node strength equal to x in the initial network Go; then after k steps of 
the algorithm, all nodes strengths will be rescaled by a factor ^ . . . , where the non-negative integers ki do satisfy 
ki + ■ ■ ■ + kg = k. Because this can be done in k\/{ki\ . . . kgl) possible different ways, we get the following relation for 
the node strength distribution for the network Gk- 

.gfc(/^.../^a:) = ^-^-^30(2;) withfci + ... + fc, = fc. (15) 

After sufficiently many steps and assuming that the main contribution arises from the choice ki ~ ■ • • fc,, ~ fc/s, we 
can use Stirling formula to get the approximate distribution (sec Fig. |4]) 

log5fc(a;) - - — ^ ^ logx , (16) 
log(/i ...fs) 

so once again the nodes strength distribution follows a power law. 



V. CONCLUSIONS 



In this paper we introduced a unifying framework for complex networks sharing several properties with fractal 
sets, hereby named Weighted Fractal Networks. This theory, that generalizes to graphs the construction of IFS 
fractals, allows us to build complex networks with a prescribed topology, whose main quantities can be analytically 
predicted and have been shown to depend on the fractal dimension of the IFS fractal; for instance the networks are 
scale-free with exponent the fractal dimension. Moreover the weighted fractal networks share with IFS fractals, the 
self-similarity structure, and are explicitely computable examples of renormalizable complex networks. 

These networks exhibit the small-world property. In fact the average shortest path increases logarithmically with 
the system size ((TT|) . hence it is small as the average shortest path of a random network with the same number of 
nodes and same average degree. On the other hand the clustering coefficient is asymptotically constant ((T^ . thus 
larger than the clustering coefRcient of a random network that shrinks to zero as the system size increases. 
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The self-similarity property of the weighted fractal networks makes them suitable to model real problems involving 
generic diffusion over the network coupled with local looses of flow, here modeled via the parameter / < 1. For 
instance one can think of electrical grids or mammalian lungs, where current or air, flows through power lines or 
bronchi-bronchioles, submitted to looses of power, or air vessels' section reduction. In all these cases the induced 
topology, namely a good choice of / and s, allows any two random nodes, final power users or alveoli, to be always 
at finite weighted distance, whatever their physical distance is, and thus to be able to transport current or oxygen in 
finite time. 



APPENDIX A: COMPLEMENTARY MATERIAL 



1. Computation of A^,"''' 



Let Ofc be the attaching node of the graph Gk- Let us define A^,"''-' ~ SigCk^'iafc' i-^- ^^'^ sum of all weighted 
shortest paths to au- Then using the recursive property and the symmetry of the map Tsj,ak "we can easily obtain a 



recursive relation for A|j"''-': 



where N^-i is the number of nodes in Gk-i- This recursion can be easily solved to get for all fc > 1 

Ai-) = (,s/)-^A^) + l-/^-M^s-l)iVo + l _ _^G./)--^-l (Al) 
^ 1 1 - / s - 1 s - f s/ - I ' ^ ' 

from which we can conclude, because / < 1, that A^°'°'' exhibits the asymptotic behavior given by equation ([5]). 



2. Computation of A^ 

Starting from the definition of the sum of all weighted shortest paths ([7]), the recursive construction and its symmetry 
we can decompose the sum Afc into three terms: 

Ak^s J2 + ^(^ - 1) E + 2^ E ^ (A2) 

where the first contribution takes into account all paths starting from and arriving to nodes belonging to the same 
subgraph, that using the symmetry can be chosen to be G^^"*. The second term takes into account all the possible 
paths where the initial point and the final one belong to two different subgraphs, and still using the symmetry we can 
set them to G^^^ and G^,^^ and multiply the contribution by a combinatorial factor s(s — 1). Finally the last term is 
the sum of all paths arriving to the attaching node ak ; once again the symmetry allows us to reduce the sum to only 
one subgraph, say G^^\ and multiply the contribution by 2s. 

Using the scaling mechanism for the edges, the first term in the right hand side of equation (|A2p can be easily 
identified with 



E ^ - 



By construction, each shortest path connecting two nodes belonging to two different subgraphs, must pass through 
the attaching node, hence using p[j^ = p|^^ + p^^^^ the second term of equation (jA2p can be split into two parts: 

E p^'- E^S^^^+ E^S^^ 

.6G1^',,6G1:^> .SGL" ,6G<f> 

where Nf, denotes the number of nodes in the subgraph G\ . Using the symmetry of the construction, the previous 
relation can be rewritten as 

E ^ = 2A^^^ E ■ 
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The last term of equation (jA2p can be related to A^°^-^ by observing that each path arriving at a^^ must pass 
through a^*^ for some iG{l,...,s}, thus 



(1) 



= 7vW + /AtT^', (A3) 

Observing that G\^^ has as many nodes as Gk-i we can conclude that Nj^^^ — Nk-i and finally to rewrite equation (jA2p 
as: 

Afc = s/Afc_i + 2s[{s - l)Nk-i + l][Nk-i + /A^t'^] ■ 
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We hereby present the construction for undirected networks, but it can be straightforwardly generalized to directed graphs 
as well. 

A straightforward generalization will be presented in the next Section IIVI See also [31 where the WFN theory will be 
generalized as to include a stochastic iteration process. 

Without loose of generality we can assume that for all integers i,j £ {l,...,m} and fc > we have f'^Xj ^ Xi and 

fifXj + l)ytXi. 



